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Abstract. Fluid dynamic limit to compressible Eulcr equations from com- 
pressible Navier-Stokes equations and Boltzmann equation has been an active 
topic with limited success so far. In this paper, we consider the case when the 
solution of the Euler equations is a Riemann solution consisting two rarefaction 
waves and a contact discontinuity and prove this limit for both Navier-Stokes 
equations and the Boltzmann equation when the viscosity, heat conductivity 
coefficients and the Knudsen number tend to zero respectively. In addition, 
the uniform convergence rates in terms of the above physical parameters are 
also obtained. It is noted that this is the first rigorous proof of this limit for 
a Riemann solution with superposition of three waves even though the fluid 
dynamic limit for a single wave has been proved. 
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1. Introduction. This paper concerns the fluid dynamic limit to the compressible 
Euler equations for two physical models, that is, the compressible Navier-Stokes 
equations and the Boltzmann equation. In the first part, we consider zero dissipation 
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limit of the compressible Navier- Stokes system for viscous and heat conductive fluid 
in the Lagrangian coordinates: 

v t ~ u x = 0, 

, u x 

ut+Px = £(—)x, / 1; q 
/ , u 2 , e x uu x 

[e + —) t + (pu) x = (k he ) x , 

2 1 v v 

where the functions v(t, x) > 0,u(t,x),9(t,x) > represent the specific volume, 
velocity and the absolute temperature of the gas respectively. And p — p(v, 9) is 
the pressure, e = e(v,0) is the internal energy, e > is the viscosity coefficient, 
K > is the coefficient of the heat conductivity. Here, both e and ft are taken as 
positive constants. And we consider the perfect gas where 

P= — =A«-Texp(l— a), e= -, 1.2 

v x R ' 7 — 1 

with s denoting the entropy of the gas and A, R > , 7 > 1 being the gas parame- 
ters. 

Formally, as the coefficients k and e tend to zero, the limiting system of (1.1) is 
the compressible Euler equations 

ih — u x = 0, 

u t + Px = 0, (1.3) 
(e + £)t + {pu) x = 0. 

The study of this limiting process of viscous flows when the viscosity and heat 
conductivity coefficients tend to zero, is one of the important problems in the theory 
of the compressible fluid. When the solution of the inviscid flow is smooth, the zero 
dissipation limit can be solved by classical scaling method. However, the inviscid 
compressible flow usually contains discontinuities, such as shock waves and contact 
discontinuities. Therefore, how to justify the zero dissipation limit to the Euler 
equations with basic wave patterns is a natural and difficult problem. 

Keeping in mind that the Navier-Stokes equations can be derived from the Boltz- 
mann equation through the Chapman-Enskog expansion when the Knudsen number 
is close to zero, we assume the following condition on the viscosity constant e and 
the heat conductivity coefficient k in the system (1.1), cf. also [17]: 

{« = 0{e) as e — > 0; 

v = K ^ > c > for some positive constant c, as e — > 0. 

Now we briefly review some recent results on the zero dissipation limit of the 
compressible fluid with basic wave patterns. For the hyperbolic conservation laws 
with artificial viscosity 

1H + f(u)x = eu xx , 

Goodman- Xin [9] verified the viscous limit for pieccwise smooth solutions separated 
by non-interacting shock waves using a matched asymptotic expansion method. 
For the compressible isentropic Navier-Stokes equations, Hoff-Liu [12] first proved 
the vanishing viscosity limit for piecewise constant solutions separated by non- 
interacting shocks even with initial layer. Later Xin [30] obtained the zero dissipa- 
tion limit for rarefaction waves and Wang [28] generalized the result of Goodmann- 
Xin [9] to the isentropic Navier-Stokes equations. 
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For the inviscid limit of the full compressible Navier-Stokes equations (1.1), Jiang- 
Ni-Sun [17] justified the zero dissipation limit of the system (1.1) for centered rar- 
efaction waves. Wang [29] proved the zero dissipation limit of the system (1.1) 
for piecewise smooth solutions separated by shocks using the matched asymptotic 
expansion method introduced in [9]. Recently, Xin-Zeng [31] considered the zero 
dissipation limit of the system (1.1) for single rarefaction wave with well prepared 
initial data and obtained a uniform decay rate in terms of the dissipation coeffi- 
cients. And Ma [22] obtained the zero dissipation limit of a single strong contact 
discontinuity in any fixed time interval with a decay rate. 

However, to our knowledge, so far there is no result on the zero dissipation limit 
of the system (1.1) for superposition of different types of basic wave patterns. In the 
first part of this paper, we investigate the fluid dynamic limit of the compressible 
Navier-Stokes equations when the corresponding Euler equations have the Riemann 
solution as a superposition of two rarefaction waves and a contact discontinuity. For 
this, we need to study the interaction between the rarefaction waves and contact 
discontinuity. 

In the second part of the paper, we study the hydrodynamic limit of the Boltz- 
mann equation [2] with slab symmetry 

ft + = ~Q(f,f), (Mx,£)eRxR+xRxR 3 , (1.5) 

£ 

where £ = (^1,^2,^3) € R 3 , f{t,x,£) is the density distribution function of particles 
at time t with location x and velocity £, and e > is called the Knudsen number 
which is proportional to the mean free path. Remark that the notation e here is 
same as the viscosity of the compressible Navier-Stokes equations (1.1), but it has 
different physical meanings from (1.1) in different equations and related contexts. 

For monatomic gas, the rotational invariancc of the particles leads to the following 
bilinear form for the collision operator 

Q(f,g)(0 - \ J J (/(£')<?(£) + /(£MO /KMC) - /&M0) 

where £J, are the velocities after an elastic collision of two particles with velocities 

before the collision. Here, 9 is the angle between the relative velocity £ — 
and the unit vector T in = {r g S 2 : (£ - £*) • T > 0}. The conservation 
of momentum and energy gives the following relation between the velocities before 
and after collision: 

? = £-[(£-£*) -r] r, 
£=£* + [(£-£*)-r]r. 

In this paper, we consider the Boltzmann equation for two basic models, that 
is, the hard sphere model and the hard potential including Maxwellian molecules 
under the assumption of angular cut-off. For this, we assume that the collision 
kernel B(\£ — 0) takes one of the following two forms, 

B(ie-e»|,0) = |(£-k.r)| = |e-&l«»£ 

and 

B(\Z-&\,9) = \Z-t;*\^b(0), &(0)eZ 1 ([O > 7r]) > n>5. 

Here, n is the index in the potential of inverse power law which is proportional to 
r 1_n with r being the distance between two concerned particles. 
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Formally, when the Knudsen number e tends to zero, the limit of the Boltzmann 
equation (1.5) is the classical system of Euler equations 



p t + {pui) x = 0, 

(pui) t + (pu\ + p) x = 0, 

(pui) t + (puiiti)a = °> « = 2, 3, 



(1.6) 



where 



p(t,x) = 

fWi(t,x) 

u, 

p(E 



-)}t + [pui(E 



-) + pu 



l\x 



R 3 „ 



MS)f(t,x,0d£, 



R 1 



¥>i(0/(*,s,£R, i = l,2,3, 



(1.7) 



-)(t,a;) 



¥> 4 (0/(*>*>0#- 



2 JR3 

Here, p is the density, u = (ui,it2,M3) is the macroscopic velocity, E is the internal 
energy of the gas, and p = RpQ with R being the gas constant is the pressure. 
Note that the temperature is related to the internal energy by E = and 
<Pi(£)(i = 0, 1, 2, 3, 4) are the collision invariants given by 

" MO = i, 

<fi(£)=ti for 1 = 1,2,3, 

MO = ^l 2 > 

that satisfy 

/ <pi(£)Q(h,g)dt = 0, for 1 = 0,1,2,3,4. 

JR3 

How to justify the above limit, that is, the Euler equation (1.6) from Boltzmann 
equation (1.5) when Knudsen number e tends to zero is an open problem going 
way back to the time of Maxwell. For this, Hilbert introduced the famous Hilbert 
expansion to show formally that the first order approximation of the Boltzmann 
equation gives the Euler equations. On the other hand, it is important to verify 
this limit process rigorously in mathematics. For the case when the Euler equation 
has smooth solutions, the zero Knudsen number limit of the Boltzmann equation 
has been studied even in the case with an initial layer, cf. Ukai-Asano [26], Caflish 
[3], Lachowicz [18] and Nishida [24] etc. However, as is well-known, solutions of 
the Euler equations (1.6) in general develop singularities, such as shock waves and 
contact discontinuities. Therefore, how to verify the fluid limit from Boltzmann 
equation to the Euler equations with basic wave patterns becomes an natural prob- 
lem. In this direction, Yu [32] showed that when the solution of the Euler equations 
(1.6) contains only non-interacting shocks, there exists a sequence of solutions to 
the Boltzmann equation that converge to a local Maxwellian defined by the solu- 
tion of the Euler equations (1.6) uniformly away from the shock in any fixed time 
interval. In this work, the inner and outer expansions developed by Goodman-Xin 
[9] for conservation laws and the Hilbert expansion were skillfully and cleverly used. 
Recently, Huang- Wang- Yang [15] proved the fluid dynamic limit of the Boltzmann 
equation to the Euler equations for a single contact discontinuity where the uniform 
decay rate was also obtained. And Xin-Zcng [31] proved the fluid dynamic limit 
of the compressible Navier-Stokes equations and Boltzmann equation to the Euler 
equations with non-interacting rarefaction waves. About the detailed introductions 
of the Boltzmann equation and its hydrodynamic limit, see the books [4], [7] etc. 
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In this paper, we will study the hydrodynamic limit of the Boltzmann equation 
when the corresponding Eulcr equations have a Riemann solution as a superposi- 
tion of two rarefaction waves and a contact discontinuity. More precisely, given a 
Riemann solution of the Eulcr equations (1.6) with superposition of two rarefaction 
waves and a contact discontinuity, we will show that there exists a family of solu- 
tions to the Boltzmann equation that converge to a local Maxwcllian defined by the 
Euler solution uniformly away from the contact discontinuity for strictly positive 
time as e — > 0. Moreover, a uniform convergence rate in e is also given. 

As mentioned above for the compressible Navier-Stokes equations, we also need 
to study the detailed wave interactions through this limiting process. 

For later use, we now briefly present the micro-macro decomposition around the 
local Maxwcllian defined by the solution to the Boltzmann equation, cf. [19] and 
[21]. For a solution f(t,x,£) of the Boltzmann equation (1.5), set 

f(t,x,£)=M(t,x,0 + G(t,x,0, 

where the local Maxwcllian M(4, x, £) = Mu^i (£) represents the macroscopic 
(fluid) component of the solution, which is naturally defined by the five conserved 
quantities, i.e., the mass density p(t,x), the momentum pu(t,x), and the total 
energy p(E + ^\u\ 2 )(t, x) in (1.7), through 

M = M[ p u 9 \ it, x, fl = — -? ( *' X) e" , (i.g) 

And G(t, x, £) being the difference between the solution and the above local 
Maxwcllian represents the microscopic (non-fluid) component. 

For convenience, we denote the inner product of h and g in L^(R 3 ) with respect 
to a given Maxwellian M by: 

JR3 M 

If M is the local Maxwellian M defined in (1.8), with respect to the correspond- 
ing inner product, the macroscopic space is spanned by the following five pairwise 
orthogonal base 

»(«) s in 

{Xi,Xj) = Sij, i,j = 0,1,2,3,4. 

In the following, if M is the local Maxwellian M, we just use the simplified nota- 
tion (•, •) to denote the inner product (•, -)m.- The macroscopic projection P and 
microscopic projection Pi can be defined as follows 

4 

Pah = ^2{h,Xj)Xj, 

3=0 

Pi/i = h- P a h. 
The projections Po and Pi arc orthogonal and satisfy 

P P =P ,PiPi =Pi,P Pi =PiPo =0. 
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Note that a function h(£) is called microscopic or non-fluid if 

MOvi(0« = 0, i = 0,l,2,3,4, 



where </?j(£) is the collision invariants. 

Under the above micro-macro decomposition, the solution f(t, x, £) of the Boltz- 
mann equation (1.5) satisfies 

P / = M, Pi/ = G, 
and the Boltzmann equation (1.5) becomes 

(M + G) t + fc(M + G) x = -[2Q(M,G) + Q(G,G)]. (1.9) 

£ 

By multiplying the equation (1.9) by the collision invariants (fi(£)(i = 0, 1,2,3,4) 
and integrating the resulting equations with respect to £ over R 3 , one has the 
following fluid- type system for the fluid components: 

p t + {pui) x = 0, 

(p Ul ) t + (puj +p) x = - J gG x d£, 

(pui) t + {puiUi) x = — J £i£iG x d£, i = 2, 3, C 1 - 10 ) 
[P(E + ^)]t + \p*x{E + M!) +pu 1 ) x = J ^iKfGxde- 

Note that the above fluid-type system is not closed and one more equation for the 
non-fluid component G is needed and it can be obtained by applying the projection 
operator Pi to the equation (1.9): 

G t + PifeM*) + Pi&G*) = - [L M G + Q(G, G)] . (1.11) 

e 

Here Lm is the linearized collision operator of Q(f, /) with respect to the local 
Maxwcllian M: 

L M h = 2Q(M, h) = Q(M, h) + Q(h, M). 
Note that the null space 91 of Lm is spanned by the macroscopic variables: 

Xj (£), .7=0,1,2,3,4. 
Furthermore, there exists a positive constant <Jo > such that for any function 
h(0 G cf. [10], 

< h,L m h >< -cr < v(\£\)h,h >, 
where f(|£|) is the collision frequency. For the hard sphere model and the hard po- 
tential including Maxwcllian molecules with angular cut-off, the collision frequency 
^(|£|) has the following property 

o<^o<Ki£i)<c(i + i£ir°, 

for some positive constants i^o, c and < Ko < 1. 

Consequently, the linearized collision operator Lm is a dissipative operator on 
L 2 (R 3 ), and its inverse L M * exists in 

It follows from (1.11) that 

G = eL M 1 [P 1 (6M :c )]+n, (1.12) 

with 

n = L M 1 [e(G t + Pi(fiG a )) - Q(G, G)]. (1.13) 
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Plugging the equation (1.12) into (1.10) gives 
p t + (pui) x = 0, 

{pu\)t + {pu\+p) x = —{p{9)ui x ) x - I ^IL^, 



(pui) t + {puiUi) x = e{p{9)u ix ) x 



£i&n x d£, i = 2,3, 



[P(E + ^-)]t + [^(E + lA-) +P u 1 ] x = e(X(6)9 x ) x + -±(p(9) Ul u lx ) x 

3 



+6 



i=2 



(p(9)uiU ix ) x - J ^i\£,\ 2 U x d^, 



(1.14) 

where the viscosity coefficient p(9) > and the heat conductivity coefficient A((9) > 
are smooth functions of the temperature 9. Here, we normalize the gas constant 
R to be | so that E = 9 and p = \p9. The explicit formula of p(9) and X(9) can 
be found for example in [5], we omit it here for brevity. 

Since the problem considered in this paper is one dimensional in the space variable 
x £ R, in the macroscopic level, it is more convenient to rewrite the equation (1.5) 
and the system (1.6) in the Lagrangian coordinates as in the study of conservation 
laws. That is, set the coordinate transformation: 



P(t,y)dy, 



t => t. 



We will still denote the Lagrangian coordinates by (t, x) for simplicity of notation. 
Then (1.5) and (1.6) in the Lagrangian coordinates become, respectively, 



and 



ft fx~\~ fx 

V V 



v t - ui x = 0, 
u\ t +p x = 0, 
Uu = 0, i = 2,3, 



-Q(fJ), 



(1.15) 



(1.16) 



+ JJ-) t + (p Ul ) x = 0. 



Also, (1.10)-(1.14) take the form 



v t - Ul x = 0, 

uit + Px = - J £,lG x d£, 

u it = £i&G x d£, i = 2,3, 

\u\i , . n 



+ L ^-) t + (pui) a 



G t - — G x + -PiKiMx) + -Pi(CiGx) - -(L M G + Q(G, G)), 

v v v e 



with 



G = eL M 1 (ip 1 (aM ;c ))+n 1 . 
Ill = L^[e(G t - —G x + iPi(6G a )) - Q(G, G)] 



(1.17) 

(1.18) 

(1.19) 
(1.20) 
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and 



V t - til 

u 



0, 



lt+Px = — Mk)x - / ZlTlixd^ 
3 v J 



I'll 



fi^ntede. * = 2,3, 



+ ^-) t +Ws = £( 



AW, 



+ey"(— 



i=2 



4e «(0) 
3 f 

:6 i' J H:,<. 



(1.21) 



With the above preparation, the main results in this paper for both the com- 
pressible Navier-Stokes equations and the Boltzmann equation will be given in the 
next section. And the proof of the zero dissipation limit for the compressible Navier- 
Stokes equations will be given in Section 3 while the proof of hydrodynamic limit 
for the Boltzmann equation will be given in the last section. 



2. Main results. 

2.1. Compressible Navier-Stokes equations. It is well known that for the Eu- 
ler equations, there are three basic wave patterns, shock, rarefaction wave and 
contact discontinuity. And the Riemann solution to the Euler equations has a basic 
wave pattern consisting the superposition of these three waves with the contact 
discontinuity in the middle. For later use, let us firstly recall the wave curves for 
the two types of basic waves studied in this paper. 

Given the right end state (v+, it+, #+), the following wave curves in the phase 
space (v, u, 9) are defined with v > and 9 > for the Euler equations. 

• Contact discontinuity wave curve: 



CD(v + , u + , 6+) = {(v,u,t 
• i-Rarefaction wave curve (i = 1,3): 



u+,P = P+,v^v + }. 



(2.1) 



Ri(v + ,u + ,9 + ) := < (v,u,6) 



v < v + , u = u + 



Xi(rj,s + )dr/, s(v,i 



(2.2) 

where s+ = s(v+,9 + ) and X{ = Xi(v,s) is i-th characteristic speed of the Euler 
system (1.3) or (1.16). 

Accordingly, when we study the Navier-Stokes equations, the corresponding wave 
profiles can be defined approximately as follows, cf. [16], [30]. 

2.1.1. Contact discontinuity. If (v— , U— , 6— ) G C'D(v+,u+,9 + ), i.e., 

u- = u + , p- =p+, 

then the following Riemann problem of the Euler system (1.3) with Riemann initial 
data 

(v-,U-,6-), x<0, 
(v+,u+,6+), x>0 
admits a single contact discontinuity solution 



(v,u,6)(t = 0,x) = 



(v cd ,u cd ,9 cd )(t,x) 



(v-,u-,6-), 
(v+,u + ,9 + ), 



x < 0, t > 0, 
x > 0, t> 0. 



(2.3) 
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As in [14], the viscous version of the above contact discontinuity, called viscous 
contact wave (V CD , U CD , Q CD )(t, x), can be defined as follows. Since we expect 
that 

P CD ^p+=P-, and |U CD |«1, 
the leading order of the energy equation (1.1)3 is 

-6 t +p+U x = k(^-) x . 



7 _1 * ■ " ■* "V 
Thus, we can get the following nonlinear diffusion equation 

®t = ae( — ) Xl B(t,±) = 0±, a= — , 

which has a unique self-similar solution QCt, x) = Q(rj), i] = . x 

V£(l+t) 

Now the viscous contact wave (V CD , U CD 7 Q CD )(t 7 x) can be defined by 

P+ 

U CD {t,x)=u + + K{l - 1)@ ^ X \ (2.4) 
Rl Q(t,x) 

fc> (t, x) = B(r, x) H fc)f. 

7P+ 

Here, it is straightforward to check that the viscous contact wave defined in (2.4) 
satisfies 

\Q CD -0 ± \ + [e(l+tp\e° D \+e(l + t)\e°°\ = 0(l)S CD e-^, (2.5) 

as \x\ — > +00, where S CD = \6 + — #_| represents the strength of the viscous contact 
wave and Co is a positive generic constant. Note that in the above definition, the 
higher order term £ ^' R7 ~^ 7 ~ 1 ^ 6 t is used in Q CD {t, x) so that the viscous contact 

wave (V , U CD 7 Q CD )(t,x) satisfies the momentum equation exactly. Precisely, 
(V CD ,U CD ,Q CD )(t,x) satisfies the system 

( v t CD -u x JD = o, 

ttCD 

jjCD , pCD _ e .(^x_\ 

"t + F x - e (T^CD>^ (2.6) 
n r\CD . tjCDttCD .t^x \ , A u x I 



7 -l f 1 x '"^yCD" • yCD 1 "b 

RQ CD 

where P CD = — j^tt- and the error term Q CD has the property that 
V 

Q CD =0(l)S CD e(l + t)- 2 e-^J, as |z| -> +00. (2-7) 

Remark 1. The viscous contact wave (V CD , U CD , <d CD )(t, x) defined in (2.4) is 
different from the one used in [14] and [16]. Here, this ansatz is chosen such that 
the mass equation and the momentum equation are satisfied exactly while the error 
term occurs only in the energy equation. However, note that the approximate 
energy equation that the viscous contact wave satisfies is not in the conservative 
form. 
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2.1.2. Rarefaction waves. We now turn to the rarefaction waves. Since there is no 
exact rarefaction wave profile for either the Navicr-Stokes equations or the Boltz- 
mann equation, the following approximate rarefaction wave profile satisfying the 
Euler equations was motivated by [23] and [30] . For the completeness of the pre- 
sentation, we include its definition and the properties in this subsection. 

If (v—,U-,8—) G Ri(v+, u+, 0+){i = 1,3), then there exists a i-rarefaction wave 
(v r \ u r \ 9 Ti )(x/t) which is a global solution of the following Ricmann problem 

vt - u x = 0, 
u t +p x {v,0) = 0, 

J*-0 t +p(v,0)u x =0, (2.8) 

7-1 

(V U 0)(t = O X) = \ («-. U -^-).^<°: 
[V,U,V){t U,X) | ( u+ , w+ ,0 + ) ;X> O. 

Consider the following inviscid Burgers equation with Ricmann data 

Wt + ww x — 0, 

If W— < w + , then the above Riemann problem admits a rarefaction wave solution 

w -, f < w~, 

w r (t,x) = w r {-) = { f, w-<f<w + , (2.10) 



t 



t > — t 

w+, f > w+. 



Obviously, we have the following Lemma, 

Lemma 2.1. For any shift to > in the time variable, we have 

C 

\w r {t + t ,x) -w r (t,x)\ < -t , 
where C is a positive constant depending only on w± . 

Remark that Lemma 2.1 plays an important role in the wave interaction estimates 
for the rarefaction waves. 

As in [30], the approximate rarefaction wave (V R , U R , & R )(t, x) to the problem 
(1.1) can be constructed by the solution of the Burgers equation 

( w t + ww x = 0, 

< ,„ s , s ,x. w + +W- w + - w_ x (2.11) 
| w(0,x) =w a {x) =w{-) = — y — + — tanh-, v ; 

where a > is a small parameter to be determined. Note that the solution w^(t, x) 
of the problem (2.11) is given by 

w^(t,x) = w a (x (t,x)), x = x Q (t,x) + w (T {x Q (t,x))t. 

And w r a (t,x) has the following properties: 

Lemma 2.2. (|30]j Let w_ < w + , (2.11) has a unique smooth solution w^.(t,x) 
satisfying 

(1) W-< w r a (t,x) < w+, (w r a ) x {t,x) > 0; 

(2) For any p (1 < p < +oo), there exists a constant C such that 
||JL<(i,.)|j LP(R) <Crnin{( U , + -^)a- 1+1 /p, (w+ - W _) 1/p ^ 1+1/p }, 
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)-t < and W- > 0, the 



\w r a {t,x) — u>-\ < (w+ — ut-)e 



21* — ui_t| 



21*— u,_t| 



2|x-ra , tl 



If x — w+t > and ui+ < 0, then 

\w r a (t,x) — w+\ < (w+ — W-)e~ 
|£<(t,.)|< 2fe^l e -^ , 
#J sup - w r (f )| < min {(«;+ - «;_), f [ln(l + 1) + | ln<r|]}. 

Then the smooth approximate rarefaction wave profile denoted by 
(yR i) U Ri ,e Ri ){t, x) (i = 1, 3) can be defined by 



s H > (t, x) = s(y H > (t, x ), e Hi (t, x)) 

w± = X l± := Xi{v±,9±), 
w'(t + to,x) = )*(y R, (t,x),s + ), 

U R ' (t, x) = u+ - 



V Ri (t,x) 



(2.12) 



Xi(v, s + )dv, 



v+ 



where to is the shift used to control the interaction between waves in different 
families with the property that to — > as e — >• 0. In the following, we choose 



to = £ 5 , and cr = s 5 . 
Note that (V R %U Ri ,Q Ri )(t,x) defined above satisfies 

?' - us* 

jRi 



0, 

< Uf* + P Ri = 0, 



(2.13) 



(2.14) 



where P Ri =p{V Ri ,@ Ri ). 

By Lemmas 2.1 and 2.2, the properties on the rarefaction waves can be summa- 
rized as follows. 

Lemma 2.3. The approximate rarefaction waves (V Ri , U Ri , Q Ri )(t, x) (i = 1,3) 
constructed in (2.12) have the following properties: 

(1) U R * (t, x) > for x £ R, t > 0; 

(2) For any 1 < p < +oo, the following estimates holds, 

\\(v R %u R %e R %\\ LP{dx) < c(t + t Q y 1+ K 

\\(V R >,U R ',Q R >) xx \\ LP{dx) < Ca' 1+ Ht + to)" 1 , 
\\{V R %U R ^Q R % xx \\ Lndx) < Ca- 2+ Ht + t r\ 

where the positive constant C only depends on p and the wave strength; 

(3) Ifx > X 1+ (t + t ), then 



| (V Rl , U Rl , 6 Kl ) (t , x) - (v-,U-,6-)\ < Ce 

\(v R ^u R ^e R %(t, x )\<^e- 2}x -^ lt+t0)} 



2|x-A 1+ (f + t ) 
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If x < X 3 -(t + t ), then 

\(v R s,u R s,e R s)(t,x) - (v+,u+,e + )\ < Ce - 2 ^- (t+t °\ 
|(y^,^,e^).(t, a )|<g e - "-^ (>HD) ' ; 

(4) There exist positive constants C and ao such that for a € (0,<7o) and t, to > 0, 

sup \(V R *,U R >,e R ')(t,x) - (v r %u r ',6 r ')(j)\ < f [<7ln(l + 1 + f„) + <r| lntr| + i ]. 

2.1.3. Superposition of rarefaction waves and contact discontinuity. In this sub- 
section, we will define the solution profile that consists of the superposition of 
two rarefaction waves and a contact discontinuity. Let (v_,u_,#_) £ R\-CD- 
Rs{v+, u+, #+). Then there exist uniquely two intermediate states w„, and 
(v*,u*,8*) such that («_ , u_ , 0_ ) g i?i(v,,u,,0 t ), (v., u», 0,) G CD(u*, u*, 6>*) and 

(«*,u*,0*) eR 3 (v + ,u+,e ± )._ 

So the wave pattern (V,U,Q)(t,x) consisting of 1-rarefaction wave, 2-contact 
discontinuity and 3-rarefaction wave that solves the corresponding Riemann prob- 
lem of the Euler system (1.3) can be defined by 



(t, x) = ( u ri + u cd + u r3 I (t, x) 

? ri + Qcd + r 3 

where (v Tl ,u Tl ,8 ri )(t,x) is the 1-rarefaction wave defined in (2. 
state (v+,u+,6+) replaced by (u*, it*, #*), (v cd ,u cd ,8 cd )(t,x) is the contact discon- 
tinuity defined in (2.3) with the states (y-,U—,8-) and (v+,u+,9+) replaced by 
(v* , u* , ) and (v*,u*,9*) respectively, and (v r3 ,u r3 7 8 r3 )(t, x) is the 3-rarefaction 
wave defined in (2.8) with the left state (v-,U-,9-) replaced by (v* ,u* ,8*). 

Correspondingly, the approximate wave pattern (V, U, 0)(i, x) of the compress- 
ible Navier-Stokes equations can be defined by 

, y Rl + yCD + yR 3 \ /«.+«* \ 

(t,x)=\ U Rl +U CD + U R ' 3 ](t,x)-\ u*+u* , (2.16) 

V e R i + e CD + e R - J V e. + e* J 

where (V Rl , U Rl , Q Rl )(t, x) is the approximate 1-rarefaction wave defined in (2.12) 
with the right state (v+,u+,8 + ) replaced by (u*, it*, 0*), {V CD , U CD , @ CD ){t, x) is 
the viscous contact wave defined in (2.4) with the states u~, 8-) and (v+, u+, 6+) 
replaced by (u» , , 0* ) and (u* , it* , 0* ) respectively, and (V R3 , C/^ 3 , R3 ) (t , x) is the 
approximate 3-rarefaction wave defined in (2.12) with the left state (i>_ , U-, 8-) re- 
placed by (v*,u*,8*). 

Thus, from the construction of the contact wave and Lemma 2.3, we have the 
following relation between the approximate wave pattern (V, U, Q)(t, x) of the com- 
pressible Navier-Stokes equations and the exact inviscid wave pattern (V, U, 0)(i, x) 
to the Euler equations 

\(V,U,Q)(t,x) - (V,U,Q)(t,x)\ 

< j[a\n(l+t + h)+a\\na\+t } + CS CD e-^, (2 ' 1?) 

with in = £ r > and a = . 
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Moreover, (V, U, Q)(t,x) satisfies the following system 

'v t -u x = o, 



U t + P x = s(^) x 



(2.18) 



R 
7-1 

where P = p(V, 8), and 



e t + pu x 



Ui 



= (P - P Kl - P CD - P H % - e(— 



)x + £ y 



u, 



u 



CD 



= (PU X - P Rl U R * - P UU U^ U - P K3 U? S ) - k(— - 



CDttCD 



V V CD 



)•<■ • 



e 



u 2 



V V CD ' 



■) 



Q 



Direct calculation shows that 

= o(i)\\(v r \q Ri )\\(v cd -v*,q cd -e*,v R3 - v *,e R3 -i 
+\(v x R3 ,® x 3 )\\(v Rl -^,e fli -e*,v CD -v*,e CD -e*)\ 



-\(v x CD ,e° D ,ug D )\\(v R > 



+0{l)e\\U. 
■= Qn + Qi2- 

Similarly, we have 

Q 2 = 



Hi I 



\Ui 



r R \V R *)\+e\{U R \- 

u R W Rl \ + \u R *\\v R3 \ 



+e\{U^ D ,V^ D )\\{U^\V^\U^,V^)\ + E\{U^\V^^,V^)\ 



(2.19) 



1\ 



O(l) \U Rl \\(V CD - v x ,Q CD - 6^,V R3 -v*,Q R3 

+\U R3 \\(V Rl -v*,Q Rl ~6*,V CD ~v*,Q CD -0*)\ 

+\(uZ D ,V x aD ,eC D )\\(V R i-v*,Q R i-6*,V R3 -v*,Q R3 

M(U^ D ,V x CD ,Q^)\\(U R \V R \Q R \U R3 ,V R3 ,e R ^\ 

+e \(u Ri , v R > , ) 1 1 (u R3 , v R3 , e R3 ) | 



-0(l)e 



le 



Hi 



^CD\ 



■■= Q 2 i + 2 2 + |Q t 

(2.20) 

Here Qn and Q21 represent the interactions coming from different wave patterns, 
Q12 and Q22 represent the error terms coming from the approximate rarefaction 
wave profiles, and Q CD is the error term defined in (2.7) due to the viscous contact 
wave. 

Firstly, we estimate the interaction terms Qn an d Q21 by dividing the whole 
domain SI = {(t,x)\(t,x) G R + x R} into three regions: 

0_ = {(t,x)\2x < \ u {t + t )}, 

Vcd = {(t,x)\\i*(t + t ) < 2x < X* 3 (t + t )}, 



{(t,x)\2x> X* 3 (t + t )}, 



where Ai* = Ai and = X 3 (v* 7 9*). 

Now from Lemma 2.3, we have the following estimates in each section: 

• In Sl_, 



698 



FEIMIN HUANG, YI WANG AND TONG YANG 



?- v *,V CD -v*)\ =0{l)S CD e- C[X M^ )] 

Ct (t + t ) 

= 0{l)e <= 
= 0{l)e 



• In flcD, 

\V Rl - 

\V Rs - v* 

• In 

\V Rl -v* 



I =0(l)e- 



= 0(l)e 

2|x|+2|A lt |(t + t ) 



_ 3/5 C(|x|+ t + t ) 



= 0(l)e-l Al *l £ 7 e 



2|x| + |A lt |(t + t ) 



2|x|+2A*(t + t ) 

G>(l)e 

,» _i/ 5 2M+A|Ct+to) 

0(l)e- A 3 £ X e ?7S : 



en 



••| = 0(l)e 
= 0(l)e 



2|x|+2|A 1 ,|(t + f ) 

,, . _i/ 5 2|x|+2lA lt |(t + t ) 
-| A 1*|E g e 2/5 ^ 

C[>S(t + t )] 2 

= 0(l)S C£, e Err+tj- 

= 0(l)e 

Ct (|x|+t + t ) 

= 0(l)e 



0(l)e 



_ 3/5 C(|x|+t + t ) 



Hence, in summary, we have 



|(Qii,Qai)| = 0(l)e 



-Ce 



-1/5 C(|x|+t + t ) 

e 7^ 



for some positive constants C. 

Now we consider the system (1.1) with the initial values 

{v,u,6){t = 0,x) = (V,U,Q)(t = 0,x). 

Introduce the following scaled variables 



t 



y 



(2.21) 



(2.22) 



(2.23) 



In the following, we will use the notations (v,u,9)(T,y) and (V, U, 0)(t, y) for the 
unknown functions and the approximate wave profiles in the scaled variables. Set 
the perturbation around the composite wave pattern {V, U, 0)(t, y) by 

(</>, V, C)(r, y) = (v - V, u - U,9 - 6)(r, y). 

Then the perturbation (</>, £)(t, y) satisfies the system 

(f>r ~ Ipy = 0, 



v 



7' 



6y Qy 

V 



(2.24) 



— Cr + iPUy - PUy) = 

,^C)(r = 0,y) = 0. 
And this system will be studied in Section 3 



2.1.4. Main result to the compressible Navier-Stokes equations. We are now ready 
to state the main result on the compressible Navier-Stokes equations as follows. 
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Theorem 2.4. Given a Riemann solution (V,U,Q)(t,x) defined in (2.15), which 
is a superposition of two rarefaction waves and a contact discontinuity for the Euler 
system (1.3), there exist small positive constants Sq and Sq such that if the contact 
wave strength S CD < 5q and the viscosity coefficient e < e$, then the compressible 
Navier-Stokes equations (1.1) with (1.2) and (1.4) admits a unique global solution 
(v e , u £ , 9 £ ){t, x) satisfying 

sup \{v £ ,u £ ,e £ )(t,x) - (v,u,e)(t,x)\ < c h e*, v/i>o, (2.25) 

where Eft = {(t,x)\t > h, ^°( +t > he a ,0 < a < |}, and the positive constant Ch 
depends only on h but is independent of e. 

Remark 2. Theorem 2.4 shows that, away from the initial time t = and the 

2 

contact discontinuity located at x = with the expansion rate 7(f+t)> f° r the 
viscosity coefficient e < en, there exists a unique global solution (v £ , u £ , 9 £ )(t, x) of 
the compressible Navier-Stokes equations (1.1) which tends to the Riemann solution 
(V, U, 0)(t, x) consisting of two rarefaction waves and a contact discontinuity when 
e — > and k = 0(e) — > 0. Moreover, a uniform convergence rate e~s holds on the 
set E/, for any h > 0. 

Remark 3. Theorem 2.4 holds uniformly when (t,x) € Eft for any fixed h > 
if the contact wave strength S CD and the viscosity coefficient e are suitably small. 
However, if we restrict the problem to a set Eft n {t < T} for any fixed T > 0, then 
we do not need to impose the smallness condition on the contact wave strength 
5 CD because one can apply the Gronwall inequality to get an estimate depending 
on time T rather than the uniform estimate in time. 

2.2. Boltzmann equation. We now turn to the Boltzmann equation. Similarly, 
we also define individual wave pattern, and then the superposition and finally state 
the main result in this subsection. 



2.2.1. Contact discontinuity. We first recall the construction of the contact wave 
(V CD ,U CD ,G CD )(t,x) for the Boltzmann equation in [16]. Consider the Euler 
system (1.16) with a Riemann initial data 

(v,u,6)(t = 0,x) = { (v-,u-,6 ) x<0, (2 26) 

where u± = (ui±,0,0) and v± > 1 6± > 0,iti± are given constants. It is known 
(cf. [25]) that the Riemann problem (1.16), (2.26) admits a contact discontinuity 
solution 

(v cd ,u cd ,& cd )(t,x) = { ("->«-> M> !C <0> (2.27) 
\ (v+,u + ,6 + ), x > 0, y J 

provided that 

26L 29+ 
u 1+ =m-, p- := - — =P+:=t. — ■ (2-28) 
3w„ 3v + 

Motivated by (2.27) and (2.28), we expect that for the contact wave 
(V CD ,U CD ,0 CD )(t,x), 

pCD = §c»*P + , \U CD \ 2 «1. 



TOO 



FEIMIN HUANG, YI WANG AND TONG YANG 



Then the leading order of the energy equation (1.21)4 is 

6 t +p + u lx =e(^^) x . (2.29) 
v 

By using the mass equation (1.21)i and v w we obtain the following nonlinear 
diffusion equation 



e t = e{a{6)e x ) x , a(0) = 



9p+A(fl) 
106* 



(2.30) 



From [1] and [6], we know that the nonlinear diffusion equation (2.30) admits a 



unique self-similar solution €>(?/), r/ 



tions 



V £ ( 1 + t ) 



with the following boundary condi- 



e(-oc,i) = #_, e(+oo,t) = e+. 

Let S = \6+ — 6—\. Q(t,x) has the property 

0(1)6 CD 

Q x (t,x) = e ^+^1, as z ^ ±00, 

with some positive constant c depending only on 8± . 

Now the contact wave (V CD , U CD ,G CD )(t, x) can be defined by 



(2.31) 



V 



CD 



3^' 
3p+ 



U? D =u 1+ 



2ea(6 ) t 
3p 4 



X) [/f D = 0,(z = 2,3), 



(2.32) 



ea D = e + ^ et[ 4 M( @ ) _3 A(e)] 



Note that the contact wave (V CD , U CD , Q CD )(t, x) satisfies the following system 



u 



CD 
It 



u. 



CD 



pCD = -/ 

/x(e c ^) 



4e ,/i(e c ^) 



ug D ) x + Q\ 



CD 



yCD 



as; y° D 

)xj i = 2, 3, 



(-) 



CDttCD 



(2.33) 



3 y CD 



CD 



where 



c/Jsr)^ = o(i)5°^e5(i + t)-% e -i(i+t) 



(2.34) 



4e/i(e CD 



CD 



= 0(l)5 CD e{l + t)- 2 e -^j, 
with some positive constant c > depending only on 9± . 



(2.35) 
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Remark 4. The viscous contact wave (V CD ,U CD ,Q CD )(t, x) for the Boltzmann 
equation (1.5) defined in (2.32) is different from the one used in [16]. Here, this 
ansatz is chosen such that the momentum equation is satisfied with a higher order 
error term. This is also different from the compressible Navier-Stokes equations 
where the ansatz satisfies the momentum equation exactly. But similar to the 
compressible Navier-Stokes cases, the approximate energy equation that the viscous 
contact wave satisfies is not in the conservative form. 

From (2.31), we have 

|e-0_| 

Therefore, 

\{V CD ,U CD ,e CD )(t,x) - {v cd ,u cd ,e cd ){t,x)\ = 0(l)5 CD e-^JTtj. (2.37) 

2.2.2. Rarefaction waves. The construction of the i-rarefaction wave 

(V Ri ,U Ri ,Q Ri )(t,x) (i = 1,3) to the Boltzmann equation is almost same as the 
one defined in (2.14) for the compressible Navier-Stokes equations in the previous 
section. By setting U Ri — for i = 1,3 and j = 2,3, all the properties of the 

approximate rarefaction waves (V Ri , U Ri , @ Ri )(t, x) (i — 1,3) given in Lemma 2.3 
will also be used later. 

2.2.3. Superposition of rarefaction waves and contact discontinuity. We now con- 
sider the superposition of two rarefaction waves and a contact discontinuity. Set 
(v-,u~,6-) G R\-CD-Rs{v-\-,u+, 0+). Then there exist uniquely two interme- 
diate states (w*, «*,#*) and (v*,u*,9*) such that (i>*,it*,0*) e Ri(v-,u~,9-), 
(u*,it*,0*) 6 CD(v*,u*,0*) amd («*,«*, 0*) e R 3 {v+, u + , 9+). 

So the wave pattern (V, U, 0) (t, x) consisting of 1-rarefaction wave, 2-contact 
discontinuity and 3-rarefaction wave as a Riemann solution to the Euler system 
(1.16) can be defined by 

V \ / v ri + v cd + v Ts \ /«»+«* \ 

Ui (t, x) = ul 1 + uf + u\ 3 (t, x)-\ uu+ut , 
6 / V 6 Tl + cd + 9 rs J \ 9* +6* j ( 2 -38) 

Ui = 0,(i = 2,3). 

where (v ri , w^ 1 , 9 Tl )(t, x) is the approximate 1-rarcfaction wave defined in (2.8) with 
the right state (v+,u + ,9+) replaced by (u*, iti*, 9*), (v cd ,ul d ,9 cd )(t,x) is the con- 
tact discontinuity defined in (2.27) with the states (v_,it_,#_) and (v + ,u + ,9 + ) 
replaced by (v* , u*, #*) and (v*,u*,9*) respectively, and (v V3 , u[ 3 , 9 ra )(t, x) is the 
3-rarefaction wave defined in (2.8) with the left state (v-,U-,9-) replaced by 
(v*,ut,9*). 

Correspondingly, the approximate superposition wave (V, U, 0)(i, x) can be de- 
fined by 

V \ ( V R i + V CD + V R * \ ( v. + v* \ 

Ui (t,x)= U R '+U? D + U R3 (t,x)- uu+ul , 

e J \ e Bi + e CD + e R * J \ 0* + 0* ) (2.39) 



= 0(1)S 



CD, 



2e < 1 +*) , if x < 0, 



= 0(l)6 CD e~^ T ^> 1 if x > 0. 



(2.36) 



^ = 0,(4 = 2,3). 
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where (V Rl , U Rl , Q Rl )(i, x) is the 1-rarefaction wave defined in (2.12) with the 
right state (v + ,u + ,0 + ) replaced by («*, Mi*, 0*), (V CD , U^ D , Q CD )(t, x) is the vis- 
cous contact wave defined in (2.32) with the states (u_,w_,0_) and (u+,it+,0+) 
replaced by (w*,ii*,0») and (v*,u*,6*) respectively, and {V Rs , U^ 3 , e Rs )(t, x) is 
the approximate 3-rarefaction wave defined in (2.12) with the left state (v-, it_, 6L) 
replaced by (v*,ul,9*). 

Thus, from the construction of the contact wave and Lemma 2.3, we have the fol- 
lowing relation between the approximate wave pattern (V, U, S)(t, x) of the Boltz- 
mann equation and the exact inviscid wave pattern (V, U, 0)(£, x) to the Eulcr 
equations 



\(V,U,Q)(t,x) - (V,U,Q)(t,x)\ 
C 

< — [crln(l +t + t Q ) +cr|lncr| 



cs 



CD, 



(2.40) 



with to = £ 5 and a = e s . 
Then we have 



@t + PU 
where P = p(V, 9) and 



v t - U lx = 0, 

e(»)x, i = 2,3, 
Aje lfa)iB + e 



(2.41) 



.a = e( : 



Q 2 , 



Qi = (P - P Ul - P CD - P U3 ) X - e( 



, f i(e)u lx n{Q CD )ug D 



{PU lx -P R >U Rl 



pCDjjOD 



V 



■)x-Q\ 



CD 



yCD 

A(6)8, \{Q CD )Q C X D 
^ u ix ) £ ( 77 T^cS 



V 



V 



yCD 



Direct computation yields 

Qi -- 



o(i) I (v Ri , e^ 1 ) 1 1 (v CD - v* , © CD - Q* , v R3 - v* , e % - e* 
+ 1 (v? 3 , e R * )\\(v Ri - v*,@ Rl - 6*,V CD - v*, e CD -9*)\ 
+\(v x CD , Q C X D , ug D ) 1 1 (v R > - «. , e fli - 0* , f* 3 - v* , e fla - 



+0(l)e 
Qu + Q 



\U%\ + \U*>i 

\CD\ 



P?\\V? l \ + \V?\\V**\ +\Q C 1 



CD I 



12 



(2.42) 



and 

Q-2 



= 0(1) |^ i ||(y Ci3 -^ ) e Ci3 -^ ) y ii3 -u*,e ii3 -0*)| 



CD 



+|(Etf»^e™)||(v*-«. > e* 



CD 



fl 3 



e*)\ 



v*,& R * 



+e\{U^ D ,V x CD ,Q c x D )\\{U R \V x R \Q R \U R \V x R \Q^)\ 

+s \(u Ri , v Ri , e^ 1 ) | \(u Rs , ^ > qr 3 } | 

+o(i) e [|e* I + |e* | + 1 (u R 1 , , ef , ^ , v** , ef 

Q21 + Q22 + IQ? 15 !. 



(2.43) 
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Here, Qn and Q21 represent the interaction of waves in different families, Q12 
and Q22 represent the error terms coming from the approximate rarefaction wave 
profiles, and Qf D (i = 1,2) are the error terms defined in (2.34) and (2.35) due to 
the viscous contact wave. 

Similar to the compressible Navier-Stokes equations case, for the interaction 
terms, we have 

|(Qii,Q2i)HO(l) e - C£ e 37^, (2.44) 

for some positive constants C. 

We now reformulate the system by introducing a scaling for the independent 
variables. Set 

x t 

y = -, t = - 

e e 

as in the previous section for the compressible Navier-Stokes equations. We also use 
the notations (v, u, 6)(t, y), G(r, y, £), ITi(r, y, £) and (V, U, 0)(r, y) in the scaled in- 
dependent variables. Set the perturbation around the composite wave (V, U, 0)(t, y) 

by 

(0, V, C)(r, y) = (v - V, u - U,9 - 6)(r, y). 
Under this scaling, the hydrodynamic limit problem is reduced to a time asymptotic 
stability problem of the composite wave to the Boltzmann equation. Notice that 
the hydrodynamic limit proved here is global in time compared to the case on shock 
profile studied in [32] which is locally in time. 

From (1.21) and (2.42), we have the following system for the perturbation (<f>, ip, £) 

( 4>t~ Iply = 0, 

^ + (P - P)y = - 3 (^ ~ - J - eQ u 

**r = - ^^) y - / te^ydt, i = 2, 3, 

Cr + ( P U ly - PU ly) = (« - ^ + ft^k) 

3 u(9)u 2 3 r r \f\ 2 

+ / a&n ly d£- / Ci^-n^de-eUiQi-eQa, 

i=2 i=l J J 

(2.45) 

where the error terms Qi (i = 1, 2) are given in (2.42) and (2.43) respectively. 

We now derive the equation for the non-fluid component G(r, y, £) in the scaled 
independent variables. From (1.18), we have 

G r - ^G y + -Pi(£iM y ) + -P^&GJ = L M G + Q(G, G). (2.46) 

V V V 

Thus, we obtain 

^[PifeM^l+nx, (2.47) 



v 

and 

Hi(r,y t Z) =L M 1 [G T -^G, + ^P 1 (aG,)-Q(G,G)]. (2.48) 

Let 

Go(r,»,0 = ^^{PiKiC^^©* + £ • ^)M]}, (2.49) 

and 

d (r, y, £) = G(r, y, £) - G (r, y, £). (2.50) 
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Then Gi(r, y, £) satisfies 

2^T l[a ! — 20 



ui „ 1. 



+^G. y - -Pi(frG„) + Q(G, G) - G 0T . 

Notice that in (2.50) and (2.51), Go is subtracted from G because || (O^, [7^)11 
(1 + £^t) -1 / 2 is not integrable globally in r. 

Finally, from (1.15) and the scaling transformation (2.23), we have 

v v 

The estimation on the fluid and non-fluid components governed by the above 
systems will be given in the last section. 



fr-—fy + -fy = Q(f,f). (2.52) 



2.2.4. Main result to Boltzmann equation. With the above preparation, we are now 
ready to state the main result on the Boltzmann equation as follows. 

Theorem 2.5. Given a Riemann solution (V,U,Q)(t,x) defined in (2.38), which 
is a superposition of two rarefaction waves and a contact discontinuity to the Euler 
system (1.16), there exist small positive constants 6q, Sq and a global Maxwellian 
= M[„ t _ Ut t gj , such that if the contact wave strength S CD < S a , and the Knudsen 
number e < Sq, then the Boltzmann equation (1.5) admits a unique global solution 
f £ (t,x,0 satisfying 

sup \\f%t,x,0-M [ y >i}t Q ] (t,x,0\\ L 2 ( i ) <C h e*, V/i>0, (2.53) 



where E/j = {(t,x)\t > h, -J== > he a ,0 < a < h}, the norm \\ ■ \\ L 2, i ) 



is 



II VM H^|(R 3 ) a7l< ^ ^ e P os tti ve constant Ch depends only on h but is independent of 
e. 

Remark 5. Theorem 2.5 shows that, away from the initial time t = and the 

2 

contact discontinuity located at x = with the expansion rate ^f+t) ' f° r Knudsen 
number e < £0, there exists a unique global solution / e (£,£,£) of the Boltzmann 
equation (1.5) which tends to the Maxwellian Mry- jj a (t, x, £) with (V", {/, 6)(£, 2) 
being the Riemann solution to the Euler equation with the combination of two 
rarefaction waves and a contact discontinuity when e — > 0. Moreover, a uniform 
convergence rate £5 in the norm L|(-^==) holds on the set E^ for any fixed h > 0. 

Remark 6. Theorem 2.5 holds uniformly on the (t,x) £ E^ for any h > if the 
contact wave strength <5 C£> and Knudsen number e are suitably small. But if we 
restrict the problem to the set E&, R {i < T} for any fixed T > 0, then we don't 
need the smallness condition on the contact wave strength 6 by using Gronwall 
inequality to get a time dependent estimate rather than the uniform estimation in 
time. 

Notations: Throughout this paper, the positive generic constants which are inde- 
pendent of T,e are denoted by c, C or C - For function spaces, H l (R) denotes the 
l-th order Sobolev space with its norm 

1 

ii/ii« = (EiiW)'' and ii-n : =imilw 
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where L 2 (dz) means the L 2 integral over R with respect to the Lebesgue measure 
dz, and z ~ x or y. 

3. Proof of Theorem 2.4: Zero dissipation limit of Navier-Stokes equa- 
tions. We will prove Theorem 2.4 about the fluid dynamic limit for the compressible 
Navier-Stokes equations to the Riemann solution of the Euler equations in this sec- 
tion. The proof is based on the energy estimates on the perturbation in the scaled 
independent variables. In fact, to prove Theorem 2.4, it is sufficient to prove the 
following theorem. 

Theorem 3.1. There exist small positive constants Si ande\ such that if the initial 
values and the contact wave strength 5 CD satisfy 

M(t)\ t=0 + S cd <S l7 (3.1) 

and the Knudsen number e satisfies e < S\, then the problem (2.24) admits a unique 
global solution (v s ,u E ,9 e ){r,y) satisfying 

su I >\(v s ,u E ,9 s )(T,y)-(V,U,e)(T,y)\<Ce*. (3 . 2 ) 
T,y 

Here A/"(r) is defined by (3.3) below. 

We will focus on the reformulated system (2.24). Since the local existence of the 
solution to (2.24) is standard, to prove the global existence, we only need to close 
the following a priori estimate by the continuity argument 

M{t)= sup ||(</>>,C)(tV)||i<X 2 , (33) 

0<t'<t v ' ' 

where x is a small positive constant depending only on the initial values and the 
strength of the contact wave. And the proof of the above a priori estimate is given 
by the following energy estimations. 

Firstly, multiplying (2.24) 2 by ip yields 



{\^) T - (p - PW V + (^ - ^ = -eQiV 



Since p- P = i?6(^ - y) + ^ and <j) T = ip y , we get 



^-^-)1P-(P-P)i> 

v V 



'2 r ' K v V' T v 

"(- - TjWy^y - eQilp ■ 

v V 



V 



v V 



y 

(3.4) 



(3.5) 



Set 

$(z) = z- 1-lnz. (3.6) 

It is easy to check that $(1) = $'(1) = and $(z) is strictly convex around z = 1. 
Moreover, 

71 7; 1 1 PV 

[JJ0$(_)] T = R&M-) - RQ(- - -)<j> T - —L<f. (3.7) 
V V v V vV 

On the other hand, note that 

' R -e*&] T = Ad - ?)C + - (3-8) 



7-1 V e /J 7-1 9^ 7-1 v e' 7-1 06 
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and 



R 



7 



1 



d-?)C 



= (1 - I )[-(pu y PU y ) + ^L )y + _ ^) - £ Q 2 ] 

- v w v ^p n"v n e )v( v v ) + e { v v> 

~^Q2 + 



, C A _ ^, 

9 { v V 



(3.9) 



~^-i(p-%- 



v6 



-4('-^)0 



K^y / # y \ 

""HS -1 - TF> + 77] - £( ^2-7{ 

6> z t> K 9 v V 9 



u y ; 



Substituting (3.7)-(3.9) into (3.5) gives 

R 



,1 1 



e$(-)] r 

7 — 1 B 



6» v w 



r' 



6> 2 



^ 1 jl _ 

~ K v V ' 



(3.10) 



where 



Ji = - e (p-P)U y -RQM^) 
Direct calculation shows that 



PVr 

vV ' 



r e r c 2 

7-1 06 



(3.11) 



,9V 



v 



J i = pu vMtz;) + 7*(^)] - [if + K-rr). + ^[(7 - !)$(-) 



l w6 



V 



V 



PUvWzz) + ^ $ (T7)] - Ilf + "(-vT)v + ^M(7 - 1)$(t7) - H^)] 



Thus, substituting (3.12) into (3.10) gives 



■v 



R ^,0., V<* 

7 9$ n - + — + ^ 

— 1 6 v v9 



where 



^ CD [*(^) + 7*(£)] + + Ktt)j/ + eQ 2 ] 1(7 - 1)*(^) - 



u 



V 



(3.12) 



(3.13) 



,6, 



1 



"^(--^y-^( 



V 



v « y ' 9 y v v ' 



. (3 : 14) 

Here, (• • • ) y represents the conservative terms which vanishes after integrating in y 
over R. 

By the strict convexity of <$>(z) around z = 1, under the a priori assumption (3.3) 
with sufficiently small \ > 0, there exist positive constants c\ and c 2 such that, 



ci0 2 < < c 2 ^ 2 , Cl C 2 < *(|) < c 2 C 2 , 

ciW 2 + C 2 )<$(^)< C2 (^ 2 + C 2 ). 



(3.15) 
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Thus, we have 

\h\dy< I (^ + ^)dy + C(r + r )- 2 \\(^0\\ 2 
+C / 5 CD e(l + er)- 1 e-^\( ( f>,C)\ 2 dy+ f e\Q 2 \\(<j>, C)\ 2 dy. 

JR JR 

Notice that the last term e | 1 1 (0, C)| 2 on the right hand side of (3.16) can be 
estimated similarly as for the terms eQiijj and eQ2g under the a priori assumption 
(3.3). Now we estimate the terms eQiip and eQ2§ on the right hand side of (3.13). 
First. 

f e|Q 1 ||V|dy= / e(|Qu| + |Qi 3 |)Mdy. 

JR JR 

From the estimation on the interaction given in (2.21), we get 



e\Qn\\tp\dTdy 
< / II^IUs / \Qu\dxdr 



JR 

-1/5 c(t+t ) , 1 

<C/ e- Ce e-^T^IIV-ll^ll^ll^r 
</?/ H^fdr + Q, / e-c^e-^^+^UW'idT 
</3 / ||^|| 2 dr + C^e- C£ 1/5 sup||V(r)||i 

[0,r] 

MyfdT + fifmpMTif+Cpe- '- 1 '', 

[0,T] 



(3.17) 



and 

/ / e|Qi 2 ||V|dTd» 

JO JR 

<e 2 / / (|K)-UK).| 2 )I^I^2/ 

JO. JR 

< / {r + ro)~ 1 U\\ 2 Uy\\idT (3.18) 

T T 

</?/ UyfdT + Cpj (T + T Q )-i\\4,\\idT 

Jo Jo 
</? / ||^|| 2 dT + 3^r ^sup|!^(r)||3 

JO [0,t] 

</3 / ||rfdr + /3 sup ||^(r)|| 2 + 

JO [0,r] 

where To = ^ = e - ^, and /3 > is a small constant to be determined later and Cp 



is a positive constant depending on j3. 
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The term eQ2g can be estimated similarly because the only difference is about 
the error term Q CD coming from the viscous contact wave in Q2- For this, we have 



: f / \Q CD \\C\dydr 

JO JR 



IICII 



(1 + er) 2 e dy 



R 



dr 



HCIII.IlCvllMl + eT)- 



drr 



(3.19) 



</3 / ||C,!rdr + ^ e 2 sup||C||E 2 

JO [0,t] 

</3||Cyl| 2 + /3sup||C||| 2 +C p ei. 

[0,T] 



(l + ET)- 2 dT 



By substituting (3.15)-(3.19) into (3.13) and choosing (3 suitably small, we can get 

~KM y )\\ 2 + \\\/(U*\U*°)(<i>,Q\\ 2 



K<t>,M(r,-)\\ 



dr 



<C I (T + T Q )- 2 \\(^C)\\ 2 dT + Ces 

(1 + er) _1 e~T^|(0, ()\ 2 dydr. 



+C6 CD e 



JR 



Now we need to estimate ||</>y|| 2 - Let v = y, then 



(3.20) 



V 



Rewrite the equation (2.24) 2 as 



(^)r " V'r - (p - P)y - EQl = 0. 



By multiplying (3.21) by ^ and noticing that 



v v 



v V 



we get 



I v v 



V 



v v 
V, 



1 1 



- (p - p)-f - i?e,(- - -) + e 



r 



Integrating the above equality over [0,r] x R in t and y, we obtain 

(i~,y)dy- 



R 



2 v v 



,Vy n 

— {^fdydr 
, j R 2v v 



C 



||(^,C,)ll 2 + e 2 ||Qi|| 



dr + C 



(3.21) 
(3.22) 



JR 



\(V yi U v ,e y )\ 2 \(^0\ 2 dyd7 



The by using the equality 



(3.23) 



Vy _ (jXy_ Vy4> 

V 



we have 



v v V v vV ' 

C-\\U 2 - |^| 2 ) < (^) 2 < C(\^y\ 2 + |^| 2 ). 



(3.24) 
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By the estimation on Qn in (2.21) and Lemma 2.3, we have 

' f T I e 2 {\Q 11 \ 2 + \Q 12 \ 2 )dydr 
Jo, jr 

(\Qn\ 2 + e 2 \(w r s ) xx \ 2 + e 2 \(w r s ) x \ 4 )dxdt (3.25) 
+ Ce 2 (t- 2 + 5-h- 1 ) 



h\\ 2 dT <C 



< C 



< Ce 

< Cei. 



JR 

-Ce- 1 ' 5 



Moreover, we have 



\(Vy,u y ,e v )\ 2 = e 2 \(v x ,u x ,e x )\ 2 

<s 2 ^ \{V x R %U^^)\ 2 +e 2 \{V^^ D ,Q c x D f 



i=l,3 

.2/. i . \-2 



= C(t + t )- 2 + CS CD e(l + eT)- x e-^. 
Substituting (3.24)- (3.26) into (3.23) gives 



(3.26) 



ll^(r,-)l| 2 + / Uy\\ 2 dT<C\\M){r,-)\\ 2 

T J ° T 

+C \Uy,Cy)\\ 2 dr + C { (T + To)- 2 \\(cPX)\\ 2 dT + Cei 



-CS 



CD 



{l + eT)- l e-^^\(<t>,C)\ 2 dydT. 



JR 



(3.27) 



Now we estimate the higher order derivatives of (i/j, (). Multiplying (2.24) 2 by — tp yy 
and (2.24) 3 by — C, yy , and then adding the resulting equations together yield 



1 



R 



+ 2( 7 - lp 



Cy]r + ^+V^ 
V V 



{(p - P) v + + [U y (- - £)] v + eQi}^ 



1 1 



vv,, t , „ ,1 1 .. ui Ul 



!>!> 
2 



(3.28) 



{(pUy PU y ) + -^(y + [VQ V {~ ~)]y + ( ^ - ^ ) + SQ 2 }(yy. 



The right hand side of (3.28) will be estimated terms by terms as follows. From 
(3.22) and (3.26), we get 



'0 JR 

< C 



(j) - P)yll)yydydT 

\{<t>v,<v)\ + \{V y ,Qy)\m0^ y v\dydr 



CL JR 

</3 ^ \\i, yy \\ 2 dr + C p I' \UyXy)\\ 2 dr + C fj / ' (r + r )- 2 !| (0, C)\\ 2 dr 
e(l + eT) -1 e"£S£|(6 C)\ 2 dydr. 



(3.29) 



JR 



+C P 5 CD 

Similar estimate holds for the term J Q T J R (pu y — PUy)^ yy dydT. 
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Notice that 



JR 



-gl/tylpyydydT 



<C / (\<t>y\ + \V y \M y \\4> yy \dydT 

/0„ JR. 



<C / (Hy\\Hyy\\Hy\\L T + 11^11^1111^11^11)^ 

T T 

<C j \\i) yy \\i\\il)y\\^\\<j) y \\dT + Ce^ / ||V^||^ to |Mt 
o. 



(3.30) 



<P ||^ y ||^T + ^(sup||0 y || 4 + £ ) / U v \\ 2 dr 

JO [0,r] JO 

<(3 f Uyyfdr + ^( X 4 + £) f Uyfdr, 
Jo Jo 

where in the third inequality we have used the fact that ||V^||l°° < Ce^ because of 
(3.26). 

Similarly, we have 

nv^CyCyydydr 
-i T v T (3.31) 

<P / WCyyfdT + Cptf + e) [ \\C v \\ 2 dr. 



The remaining terms can be estimated directly by using (3.25) and the fact that 

[Uy{\ - = 0{l)[\{Uyy,UyVy)\\4>\ + \Uy\\<f>y\], 

[v®y{\ - y)]y = 0(1)0(9^,6^^)1101 + |6 y ||^|]. 
Hence, if we take (3 suitably small, then we obtain 

ll(^Cy)(^-)H 2 + f T W^yyXyyWdr 



<C W^y^yXyWdr + C (T + T Q )- 2 \\(fr()\\ 2 dT + Ces (3.32) 



+CS on / e(l + eT)- L e-^\((j),()\ dydT. 

JO JR 

The combination of (3.20), (3.27) and (3.32) yields that 
||(0, V,C)(t, Oil? + [ T [ll^f + 11(^)1 



dr 

<C I (r + T Q )- 2 \\(cj ) ,C)\\ 2 dT + Ce^ (3.33) 
Jo T 

+C5 CD [ [ e(l + er)- 1 e-^|(0,C)| 2 ^r. 

Jo JR 

In order to close the estimate, we only need to control the last term in (3.33), 
which comes from the viscous contact wave. For this, we will apply the following 
technique by using the heat kernel motivated by [13]. 

Lemma 3.2. Suppose that h(r,y) satisfies 

h e L°°(0,+oo;i 2 (R)), h y € L 2 (0, +oo; L 2 (R)), h T G £ 2 (0, +oo; H~ 1 (R)), 
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Then 



<C a 



:(1 + st) 1 e h 2 (T,y)dydT 



\\K0,vW+ \\h y \\ A dT+ (h T ,hgi) H - lxm dT 



(3.34) 



where 

-oo 

and a > is the constant to be determined later. 

The proof of Lemma 3.2 is similar to the one given in [13]. The only differ- 
ence here is that we need to be careful about the parameter e in the estimation. 
Therefore, we omit its proof for brevity. Based on Lemma 3.2, we can obtain 

Lemma 3.3. There exists a constant C > such that if 5 and Eq are small 
enough, then we have 



e(l + £T) _1 e Ofdydr 



JR 



<q|(^,C)(T,-)|| 2 + C7 / W&y^yXy) II 2 dr 



(3.35) 



+C / {T + T a )-^\\{4>^X)\\ 2 dT + Ce>. 
Jo 

Proof. From the equation (2.24) 2 and the fact that p — P = ^~ , we have 

i'r + (" 



,RC-P<j>, ,u v U 



)y - _ -y)v _ 



Then 



(R( - P^)y = R( v P4> {Vy + <j>y) - Vi W + - ^)y - VSQ,. 



(3.36) 



Let 



G 6 (r, </)=£(! + ex)" 1 



_ ber}~ 

e 1 + c ' r dr], 



where b is a positive constant to be determined later. Multiplying the equation 
(3.36) by G b (R(-P(f>) gives 



G b (R( - P<j>) 2 ] (R( - Pcj>) 2 

- {^bjy 



G b (R( - P<j>Y 



(Vy + 4>y) - G b v(R( - P<j>)^ T 



(3.37) 



+G b v(R(-P^ Uy 
v 



-, y -eG b v{RC-P4>)Qx. 



Note that 

-G b v(R( - P<i>)4w = -[G b v{R( - P^] T + [G b v(RC - P<j>)^] v 
+ (G b v) T (R( - Pcf))i> + GMiRC - P<f>) T , 



(3.38) 
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(R( - P4>) T 

= i?Cr ~ PA - P<Pr 



(7-1) 

-jPl/jy - P T 



U 2 U 2 



1 v V 



) y - eQ 2 



(3.39) 



By using the equality 



GbvyPM = -hG b vP^] y + 7 vP(G b ) y ^ + 1 (vP) y G b ^, 



(3.40) 



we have 

s{l+ST)- l e-¥&[{RQ-P<j>f +~jPvi> 2 ] = [G b v(R(-P<j>)Tp} T + (---) y +Q 4 , (3.41) 
where 



-(vG b ) T v[R( - *W ~ ?Y-(Pv) y G b + G b vi>P T <j> 



+ (7 - l)G b vtp 



(p-P)( Uy +^ y )-C^^) + eQ, 
lyy v V 



(3.42) 



+ [G b v(R( - P<b)} y (^ - ^JL) + ( 7 - i) v [G b wP) v {^- - ^) 

v V 



G b (R( - p<t>y 



{Vy+4> y )+eG b v(RQ- Pc^Ql 



Note that 

||G 6 (t,0|| £ - <C Q e5(i + er)-i 
Thus, integrating (3.41) over (0, r) x R gives 

e(l + er)- 1 e-££ [(R( - P<j)) 2 + ^ 2 ]dydr 

< CW^, ^ C)(t, -)|| 2 + C f ||(0,, i> y , Q(r, -)fdT 
+C [ {r + r )-i\\{cpAX){r,-)\\ 2 dT + Cei 

^-\{(j>,0\ 2 dydT. 



o Jr. 



(3.43) 



-CS 



o 

CD 



e(1+et) e 



/O JR 

In order to get the desired estimate stated in Lemma 3.3, set 

( + P4> 



7-1 

in Lemma 3.2. We only need to compute the last term on the right hand side of 
(3.34) for this given function h. From the energy equation (2.24) 3 , we have 



e, 



K = -(p - P)i>y + [P T - - P)Uy] + Kf - -f)y + ~ -f) ~ SQ 2 

5 



Thus 



/*T ^ pT f 

I {h>T,hgl) H i xH -idT = V] / / hg 2 a H t dydi 

JO i=1 JO JR 



(3.44) 
(3.45) 
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By noticing that 

||fl«(r,-)||i- <C a , 



we can 



estimate / / hg^Hidydr(i = 2, • • • ,6) directly. The estimation on 
Jo Jr 

nhg^Hidydr is more subtle. Firstly, by using the mass equation (2.24) 1; we 
t 



have 

hglHx =-(p-P)i) y hg 2 a 



v Zv 

, ( 7 -l)fe 2 0gg | lPW 2 g* | 2( 7 -l)fe 2 + 7 ^ 2 

v 2v v 

2( 1 ~l)h^ + 1 Phcf 2 , iH 2 gl n , r 2 (7-i)# , 7^ 2 



9a{ga)r 



2« 2 + ~^t Pt + [ — s — + — ^ 

5 

Now the terms Jj(i = 1, • • • ,4) can be estimated directly, cf. [13]. Here we only 
calculate the term J5. From (3.44), we have 



2w 

i=l 



Now Jg can be estimated as follows: 



I^J \4\dydr <c£ J 



<C/ ||(0,C)llLll^llll(^C)IKr 

r 

<C ||(0,C) y ||||^llll(^C)l| 2 rfr 



<Csup||(0,C)(r,.)|| 2 / \\{4>,i>,Qy\VdT 

[0,r] JO 

<c x 2 rii(0,^,c)yii 2 rfT. 



Note that the other terms J§(i = 2, ■ • • ,5) can be estimated directly, we omit the 
details for brevity. 

Therefore, by taking the constant a=^,we obtain 

. , Cpey 2 

e(l + £T) e !+" h dydr 

'o Jr 



<C||(^^,C)(t,-)I| 2 + C / ||(</)„^,C,)|| 2 rfr + C / (r + T )-*||(«/),V)|| 2 rfr 

Jo Jo 

+Ce'i +C{S CD + x ) I I e{l + eT)- 1 e- S ^\{(l>,0\ 2 dydT. 



'0 JR+ 

(3.46) 
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By taking b = Co in (3.43) and by combining the estimates (3.43) with (3.46), 
we yield the desired estimation in Lemma 3.3 if we choose suitably small positive 
constants S CD , Eo and x- d 

Now from (3.33) and Lemma 3.3, if the strength of the contact wave 6 and 
the parameter x on the a priori estimate are suitably small, we can get 



l(^,0(r,-)ll? 



\M 2 + \\Wv,<y)\ 



(It 



< C 



With this, the Gronwall inequality gives 
ll(^,C)(r,-)||?- 



l(^,C»)lli 



dr < Ce*. 



And then wc complete the proof of Theorem 3.1 by Sobolcv imbedding. 



4. Proof of Theorem 2.5: Hydrodynamic limit of Boltzmann equation. In 

the last section, we will prove the fluid dynamic limit for the Boltzmann equation to 
the Riemann solution for the Euler equations as stated in Theorem 2.5. Again, the 
proof is based on energy estimates for the Boltzmann equation (2.52) in the scaled 
independent variables. For this, it is sufficient to prove the following theorem. 

Theorem 4.1. There exist two small positive constants Si, e\, and a global 
Maxwellian M* = Ml ^^ji such that if the initial data and the strength of the 
contact wave 5 satisfy 

Af(r)\ T=Q + S CD < S u (4.1) 

and the Knudsen number e < E\, then the problem (2.52) admits a unique global 
solution / £ (t, ?/,£) satisfying 

sup||/ e (r,t/,£) -M^ejCT.^OIlLlf^) < Ce ^ '• (4.2) 
Here, jV(t) is defined by (4.5) below. 

Remark 7. If we choose the initial data for the Boltzmann equation (2.52) as 

/ £ (0,y,£) = M [v ^ e] {0,y,Z) = M [mv) , I r(o,„),e(o,v)](0. ( 4 -3) 

then 



AA(r)| T=0 = 0(1) 



\\(e y ,U y )\\ 2 + \\(V yy ,Q yy ,U yy )\\ 



0(1)£5. (4.4) 



In this case, the functional measuring the perturbation A/"(t) at t = is smaller 
than the estimate given in Theorem 4.1 that is of the order of 0(ei) because e is 
small. 



Consider the reformulated system (2.45) and (2.51). Since the local existence of 
solution to (2.45) and (2.51) is now standard, cf. [11] and [27], to prove the global 
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existence, we only need to close the following a priori estimate by the continuity 
argument: 



M{t)= sup |j(^,C)(rV)lli + 

0<T'<T 



|Gif 



-d£dy 



£ 

M=i 



\d a G| 2 
M* 



-d£dy+ J2 

l«l=2 



|^/| 2 

M* 



(4.5) 



dtdy < x 2 



where d a , d a denote the derivatives with respect to y and r respectively, and x 
is a small positive constant depending on the initial data and the strength of the 
contact wave, and M* is a global Maxwcllian to be chosen later. 
Note that the a priori assumption (4.5) implies that 



and 



IIW^OIIL <c x \ 

|Gi„| a 



(4.6) 



Q2 

TT-d^dy 



<W <C(e + x"), (4.7) 



and for |a| = 1 

| 5 « G |2 



|<9« G |2 



-d£dy 



\d a G y \* 



From (1.17) and (2.41), we have 

T - "01y = 0, 



dtdy) <C(e+x 2 )- 
(4.8) 



tfir + (p - P)» = - eQi - / CiG„de, 



3 V y 



0j 



■(^Ciy)y- / ^fiGvde, i = 2,3, 



Cr + (Pttlv - PE^lv) = -(^6 



v 3 y 



(4.9) 



Thus 



|| (</>., Vr,Cr)|| 2 <C(£ + X 2 ). 

Hence, we have 

||(«r, «r, ^r)|| 2 < C\\ (<t> T ,lpr, (r)\\ 2 + C\\ (K, U T , 9 r )|| 2 < C{e + X 2 ). 

In addition, (4.5) also implies that 

\\(v y ,u y ,e y )f < C\\{<j> y ^ y ,Qf + C\\{V y ,U y ,@ y )f < C(e + X 2 ). 

Since 

\d a f\ 2 



\\d a [p,pu,p(E+^))\\ 2 <C 



(4.10) 
(4.11) 
(4.12) 
(4.13) 
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the inequalities (4.11)-(4.13) give 

\\d a (v,u,9)\\ 2 < C\\d a Upu,p(E+ l -^)j || 2 

+ c E / \° a U <*> p( E + ¥)) i 4rf y (4 - 14) 

<c( £ + x 2 ). 

Thus, for |a| = 2, we have 

||d Q (^,C)H 2 < C(\\d a (v,u,9)\\ 2 + \\d a (V,U,e)\\ 2 ) < C(e + X 2 )- (4.15) 
Finally, from the fact that / = M + G, we can obtain for \a\ = 2, 

<C J I^JLd^dy + C\\d a (v,u,d)\\ 2 + C £ J \d a '(v,u,9)\ 4 dy (^) 

\a'\=l 

<C(e + X 2 ). 

Before proving the a priori estimate (4.5), we list some basic lemmas based on the 
celebrated H-theorem for later use. The first lemma is from [8]. 

Lemma 4.2. There exists a positive constant C such that 



M U M JM JM J M 

where M can be any Maxwellian so that the above integrals are well defined. 

Based on Lemma 4.2, the following three lemmas are taken from [20]. And the 
proofs are straightforward by using Cauchy inequality. 

Lemma 4.3. If 9/2 <9± < 9, then there exist two positive constants a = a(v,u,9; 

andrjo = r) (v, u, 9; v*, u*, 6±) such that if\v-v ir \ + \u~u+\ + \9-9+\ < rj , 
we have for h(£) G 0T- 1 , 



M* J M, 

Lemma 4.4. Under the assumptions in Lemma 1^.3, we have for each h(£) € Ot -1 ", 



M* 1 1V1 1 * - J M* 

Lemma 4.5. Under the conditions in Lemma 4-3, for any positive constants k and 
X, it holds that 

— m: — di J el - Ck ' x J — m; — c ' 

where the constant Cu.\ depends on k and A. 
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With the above preparation, we are ready to perform the energy estimation as 
follows. Firstly, similar to (3.13), we can get 



^^ 2 + i?e$(-^) + eo>(|; 



3 v 



',9V, 5_,t> ' 
y vG } 3 K V> 



-pu? v d 



V 



3 V 



_$(_) _ $( — ) 

3 y v' y e' 



A Ai{6) n{Q) 
o(— —WiyWiy 



3 V v 



V 



C y ,X(9) A(9)._ , C0y.\(0)0 y A(0)e„. , 4C/z(0) 2 /z(e) 



r 



V 



(4.17) 



where 



3 3 

#i = -£>y i' • ~(I> / fc&nn^ -5/ Cil€l 2 n ly d€)- (4.i8) 



The estimation on the macroscopic terms in (4.17) is almost same as (3.20) for the 
compressible Navier- Stokes equations so that we have 



11 



(MW + II 0(^,01 



+C(5 CI? e / / (l + £T)- x e-^-((t) 2 +( 2 )dydT+ f 
Jo Jr Jo 



(4.19) 



Now we estimate the microscopic term J J Nidydr in (4.19). For this, we 
only estimate the term T\ =: — J J ^ l J £i^yd£,dydr because other terms in 



Nidydr can be estimated similarly. 
For Ti , integration by parts with respect to y and Cauchy inequality yield 



T\ = J J i' ly J gRxd&ydT 



</3 / U ly \\ 2 dT + C / / I / gn^tfdydr. 
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By (2.49), we have 



J J | J e^idi\ 2 dydr 



+C Jo /l/^ L M^ Pl ^ G y)]^| 2rf ^ T + C / / I / H 2 i^MlQ(G,G)]d{;\ 2 dydT 

i=l 

(4.21) 

Let M* be a global Maxwcllian with its state (v+, u*, 8+) satisfying ^9 < 6+ < 6 
and \v — u*| + \u — it*| + \6 — 0*| < rjo so that Lemma 4.3 holds. Then we can obtain 



v 



(4.22) 



Similarly, 



Tf < c i - rffrfrdr. (4.23) 



Moreover, 



TS^n T f \ ^(KDILm^Pi^iG,)]! 2 ^ /■ _ Vl , h ,4 A/r ,,,, . 

JO J J iV1 [2i>*,2u*,2e*] J 



(4.24) 



From Lemma 4.2, we have 



2 i in. |2\ /• IHJ2 i in. |2 

(4.25) 



< C(s+ x*)/7/«w + ^. 
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Substituting (4.20)- (4.25) into (4.19) yields that 

||(^,C)(t,-)H 2 + [ T \\\(i'yXy)\\ 2 + \\J(U^^)(<P,0\\ 

T Jo L 

<C [ (T + r a )- 2 \\(^0\\ 2 dr + Ce^ 



(It 



+C5 CD e / / (1 + er)- 1 e-^r (0 2 + C' 2 )dydr 



JR 



(4.26) 



l«'i=i 



'0 



+C(x 2 + e)l I |»^, 



epig 



To recover the term H^H 2 in the integral / ■■■dr in (4.26), as in the previ- 

Jo 

ous subsection for the compressible Navier-Stokes equations, we firstly rewrite the 
equation (2.47)2 as 



4. M (6), , 4 r ,//(0) n(Q) s . _ /V a _ ( 4 - 27 ) 

= -3(7)^ _ 3 [( ^T " v )uiy]y + sQl + J fi n i»*' 

by using the equation of conservation of the mass (2.47) 1. 
Since 

and 

(f>yi>\T = (<Py^l)r ~ (<Pr^l)y + lpl y , 

by multiplying (4.27) by <p y , we get 



M&),2 . M P ,2 ,2/^(6) 2 2 2 



Integrating (4.28) with respect to r,y and using the Cauchy inequality yield 



ll^(r,-)lr+/ U y rdr<C\\Mr,-W + C WWyXyWdr 
T Jo Jo T 

+CS CD e [ /(l+er)- 1 e -^|(0>,C)| 2 rfydr + C / (r + r o )- 2 ||(0, C)|| 2 dr 



JO JR 

+Cei +C X Hlyy\\ 2 dT + 



I I 1 1 ^vWdytr- 

(4.29) 
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For the microscopic term J J \ J $TL\yd£\ 2 dydr, by (2.50), we have 

f j\j &R ly dtfdydrr 

[ J\J e i {^G T ) y d^dydT + £ J | J e i (^ l ^Gy)y^\ 2 dydT 
J\ J e i [^lPi(CiG y )} y d^\ 2 dydr + £ J | J g[L^Q(G,G)} v d£\ 2 dydT 



< C 



t=l 



l«|=2 



Similar estimates hold for T\ (i = 2,3). Moreover, 
T| <c£ J | J £i l m 1( 3( G ) G y )d^\ 2 dydT 

+ C fJ\J ^M{Q(^Q(G^),M v )}d4\ 2 dydT 

+c r f f^m\' d( f\^l d(dtdr 



M* J M* 



JO J 

<C{x 2 +s)£ J 



KKI)(|G 1 |" + |G.|-) 



(4.31) 



(4.30) 

Note that the inverse of the linearized operator Lj^ 1 satisfies that , for any h € Af 1 - , 

(L^ftjr = Lm(^) - 2L M 1 {Q(L M 1 / l ,M T )}, 
(L^h) y = L M \h y ) - 2L M 1 {Q(L M 1 ^,M,)}. 

Then we have 

T\ <c£ J | J e^Gyrd^dydr 

+ c fJ\J ei^{Q(L^G T ,M y )}ded y dT 

[ [ f "-^\9 a GMdydr 



+c£J \(v y ,uy,6 y )\ 2 J ^K ^r 
^ ' MI£l)|G r | 2 



(4.33) 
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ll^(r,-)H 2 + / Hy\\ 2 dT<C\\Mr,-)\\ 2 + C / \\(i; y ,Cy)\\ 2 dT 

+C5 CD e [ t {l+err 1 e~ 9 ^\{^X)\ 2 dydT + C [ (r + r )- 2 ||(^ Ofdr 



l«l=2 



+C(e + X 2 ) ^ 
Jo 



Kkl)(E K |=il^'G| 2 + |G 1 | 2 ) 



d^dydi 



M 4 

(4.34) 

We now turn to the time derivatives. To estimate \\{4> T , ">Pt> Ct)|| 2 , we need to use 
the system (4.9). By multiplying (4.9)i by 4> T , (4.9) 2 by Vir, (4.9) 3 by (i = 2, 3) 
and (4.9)4 by Cr respectively, and adding them together, after integrating with 
respect to r and y, we have 



\\(4> T ^ T ,Cr)(r,-)\\ 2 dr<C / llfa^&jfdT + Ce* 

Jo 

+ J\r + T )- 2 \\(^i',C)\\ 2 dr + C J J t4m\G v \ 2 dtdydT (4.35) 

Jo JR 

The microscopic component Gi can be estimated by using the equation (2.52). 
Multiplying (2.52) by gives 

_ -P^Gy) + Q(G, G) - G 0r j^i. 

(4.36) 

Integrating (4.36) with respect to r, £ and y and using the Cauchy inequality and 
Lemma 4.2-4.5 yield that 

<Ce^+Cj2 I W'MXWdr + C I f f V -Ml\ Gy \ 2 d^dyd Tl 
I" 1=1 

where we have used the fact that 

/ J^( T = ^^dy = J l^( T = 0,y,mdy 

<c||(e„4)(r = o,-)|| 2 <c £ i 

Next we derive the estimate on the higher order derivatives. By multiplying 
(2.46) 2 by -ipiyy, (2.46) 3 by -ipi yy (i = 2,3), (2.46) 4 by -£„„, and adding them 
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together, we obtain 

i=l i=2 

3V v lyriyVtyy / A lyVwVwy \ ly^y^yy 
4 r/M (6>) M(6) Nr7 l ~ 2 , u \{0) A(6V 

4 uffl uf6) (4-38) 

3 



+ E — u % - ( £ &- - zQWKvv + E f fi&ni„d£ 

i=2 i=l ^ 

3 

/ ei^de-i / eilCI 2 ni,dO- 

i— 1 

Integrating (4.38) with respect to r,y and £ yields 

ll(V'y,C ! ,)(T,-)l| 2 + [ T \Myy,(yy)\\ 2 dT 

T Jo T 

<c/ ||(^,^,Cy)H 2 rfr + C / (T + To)- 2 ||(0,^C)|| 2 dr + Ce^ 
Jo T Jo 

+C6 CD e I [ (l + er)- 1 e-%&\(<l>,il>,Q\ 2 dvdT 

JO JR 



+C(e? + X ) 



l«'l=i 



|o|=2' 

+G( e *+x) E f I I ^r\d a 'GfdZdydT 



(4.39) 

Again, to recover ||0 TO || 2 in the time integral in (4.39), by applying d y to (2.46)2, 
we get 



Tplyr + (P- P)yy = -|(^tfl V )w - eQly ~ J gG yy dl;. (4.40) 



Note that 



(p- p )w = -^yy + ^Cyy ~ \(jP-P)Vyy - \p yy -^(p-P) y - ^</y (4.41) 
Multiplying (4.40) by —<j> yy and using (4.41) imply 

ipi y (j) yy {T,y)dy + I I ^-4> 2 yy dydT <C j \\ (ip lyy , ( yy )fdT + Ce% 



+Cj(T + T )- 2 \\(<t>,i>,0\\ 2 dT + C(s-*+x)J W&y^yXyWdT (442) 

E [J I ^OMdydr. 

\a\=2 J0 J J 

To estimate || {4> yT , tp yT , ( yT )\\ 2 and || (</> TT , -0 rT , ^ TT )|j 2 , we use the system (4.9) 
again. By applying d y to (4.9), and multiplying the four equations of (4.9) by <f> yT , 
ipiyn i'iyr {i = 2,3), Cyr respectively, then adding them together and integrating 
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with respect to r and y give 

\\{4>yr^yrX V r)\\ 2 dT <C I || {<j) yy , tp yy , Q yy ) \\ 2 d T + Ce% 

JO 

+C J T (r + T Q )- 2 \\(^i',C)\\ 2 dr + C(e^ + X ) J II (</> v , (y)\\ 2 dr 

+c l 7 / w |Gy|2rfWr+c £ / T // w |a " G|2 ^ r - 

M=2 

(4.43) 

Similarly, we have 

Wicprr^rrXrrWdT <C / 1 1 ((f> yT , Vv , Cyr ) 1 1 <*r + Cs * 

JO 

+C J T (T + T Q y 2 \\(^i,,()\\ 2 dT + C(ei+x) £ f \\d a 'MX)\\ 2 dr 

+C [II 'fW*** + c £ I// ^^W***- 

|q|=2 

(4.44) 

A suitable linear combination of (4.39) - (4.44) gives 
\\^yXy,K)(r,-)\\ 2 + E r W^Mfdr 

^ e /7 / 4P |9aG|2dWT+c Ef// 4P |aa ' G|2 ^ r 

, , 2 Jo J * |a'| = l"' J J * 



M=2" u |a'|=l 



c( £ '+x) V ri|a Q '(^>,c)|| 2 dr + c £ 

|a<| = l J ° 



|Q'| = 

(4.45) 

To close the a priori estimate, we also need to estimate the derivatives on the 
non-fluid component G, i.e., d a G, (\a\ = 1,2). Applying d y on (2.47), we have 



Since 



we have 



G yT - C^G y ) y + {ip^M,)}, + {±Pi(6G„)}„ 
LmG s + 2Q(M. y , G) + 2Q(G y , G). 



Pl(ClMy) = 2^ Pl[a( ^"^ + * ' Uy)M] ' 



(4.46) 



li-P^iM^I < C(«J + + 2 + |^| + | %y |)|s (e) | M , 

where -B(£) is a polynomial of £. This yields that 

[j J \{\^^y)}y^mydr<lfJJ ^\G y \^ d ydr 

+C f W(i>yyXyy)\\ 2 dT + C(ei + x ) f \\{^ y Xy)\\ 2 dr + Cei 
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Thus, multiplying (4.46) by and using the Cauchy inequality and Lemmas 
4.2-4.5 yield 



Kiel) , J 



<C(e*+ X ) J o J J -j^lG^dtdydr + C^ + X ) J \\{M v ,CyWdr 
+ °Io 1 1 ^\ G yy\ 2d ^ d V dT + c [ W(^Cyy)\\ 2 dr + Cel. 

(4.47) 



Similarly, 



(T,y,Z)dSdy+ / / / -^\G T \ 2 d^dydT 



2M* v J J J M* 

<C [ [ [ ^\G yT MdydT + C^ +x ) f J J^\G 



+C(ei+ X ) ^ J J ^M\ Gy Mdydr 
+C(si+x) £ T \\d a 'M,Otdr + C [ T 



where we have used the fact that 



\\(^yrXyr)\\ 2 dT + Ce^ 

(4.48) 



2M* v ^ * j J 2vM* 

<q|( V , M ,0),(r = O,-)|| 2 
= q|(T/,[/,e),(r-0,-)|| 2 <Cei 

Finally, we estimate the highest order derivatives, that is, J ipi y 4>yydy and 
Jo" / / ^m^ 1 * d t d V dT with |a| = 2 in (4.45). To do so, it is sufficient to study 

/ / l -^ff d £ d y (M = 2 ) in view of ( 4 - 13 )- ( 4 - 16 )- For this > from ( 2 - 53 ) we have 

vf r - ui/„ + £i/„ = vQ(f, f) = i>[L M G + Q(G, G)]. 
Applying <9 Q (|a| = 2) to the above equation gives 

V(d a f) r - VL^G ~ Ul (d a f)y + tl{d a f)y 

= -d a vf T + d a Ul f y - [d a ~ a 'vd a 'f T - d a - a ' Ul d a ' f y ] (4.49) 

|c'|=l 

+ [d a (vL M G) - vL M d a G] + d a [vQ(G, G)]. 

Multiplying (4.49) by g£ = ^ + ^ yie lds 

dd a f| 2 <9"G 
(^)-^G — 

= ^{ - ^/r + 9 Q «l/ y - 2 Id^'^'fT - d^' Ul d a ' fy] 
* |Q ' I = 1 

+ [d a (vL M G) - vL M d a G] + O a [vQ(G, G)]| + dLm^G ■ ™ + (■ • • ) v . 

(4.50) 
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Hence, 



<(/L|/ (l M, l + | G ,|)l™' + '«l 



<C(e + X 2 ) / \\d a {^X)\\ 2 dr+^- 

Jo 2 16 Jo 

+C(£ + X2) [I I^ dVdT 



«|<9 Q G| 



■d^dydr 



C(si+ X ) I' \\d a '(^,0\\ 2 dr + Cei, 



and 



E 

|a'| = l 

* E 

|a'| = l 
(T 



< 



16 



\d a -<*'vd a 'f T ^l\<%dydT 
f \d a - a 'v\ J (\d a 'M T \ + \d a 'G T \) ldaMl ^ daGl d^dydr 
^^Id^dydr + C(S + 7 ) f ||3 Q (</>, ^, Ofdr + Cel 



Notice that similar estimates can be obtained for the terms d a uif y ^$- and 

E|o'|=iS a - a '«ifl a 'Argf. 
Furthermore, we have 



a a (wL M G) - vL M d a G = (<9%)L M G + 2vQ(d a M, G) 

^ |2wQ(9 Q - Q 'M,a Q 'G) +9 Q " Q 'w[L M 9 Q 'G + 2g(a Q 'M,G)]|, 

|a'|=l L 



and 



d a [vQ(G, G)] = (9 Q w)Q(G, G) + 2wQ(9 Q G, G) 

^vQ(d a ~ a 'G, d a 'G) + 2(d a - a 'v)Q(d a 'G, G)] 



l«'|=i 
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For illustration, we only estimate one of the above terms in the following because 
the other terms can be discussed similarly. 



a t T f fv\d a G\ 2 



-ml J J -mT^ 



° +^+*')/7/ ^ l)l|G ^ +|Gi|a w . 

f T f f d a M 
Now wc estimate the term / / / vLm<9 q G ■ — d^dydr in (4.50). First, note 

that Pi(<9"M) does not contain the term d a (v,u,9) for |a| = 2. Thus, we have 

V L md «G. d *M <dydT= r r /w l( yM)^ 



M * " In I I M 

lop 12 



|a'|=l 

(4.51) 

Also we can get 

f J J^G. d "M(±-±Wdydr<^r J J ^^GMdydr 



+Cr,£ / \\d a {<P^XWdT + C{e* +x) 2^ / \\d a (^iPXWdT + Ces, 
Jo \ a >\=i J ° 

(4.52) 

where the small constant 770 is defined in Lemma 4.3. The combination of (4.51) 

f T f f d a M 
and (4.52) gives the estimation on J J J i;Lmc? q G • — d^dydr. 

Thus, integrating (4.50) and using the above estimates give 
<C(ei+x)T I \\d a '&^X)\\ 2 dT + 0(^ + 6 + ^^2 [ \\d a (4>, V, Ofdr 

\ a >\ = l J ° \ a \=2 J ° 



+C^+ X )J2 f I [^\d a 'Gfdt;dydT + Cei 

\a'\ = l J ° J J * 

+C(ei+ x ) J r J j ^\G^d^dydr, 
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where we have used the fact that 



{t = 0,y,£)d£dy = / — ^ (r = 0,y,£)d£dy 



2M* v >»>^ * » j j 2M, 

< c|| (v, u, e) vv (r = o, -)|| 2 + c\\(v, u, 0) V ( T - o, 



Finally, similar to Lemma 3.3 in the previous section, we can get 

C £ 2 

e(l +eT) _1 e _ T^|(0, ip, C,)\ 2 dydr 

R 

<Cj|(^V,C)(r,0H 2 + C / \\^ y ,^ y ,Q\\ 2 dT + Cei 

Jo 

+C | r (r + r o )-f||(0,^C)|| 2 dT + C( £ 5 +Y ) | T | J " m ^ Gl]2 d^dydr 



|a'|=l' 

Note that here we need to estimate the microscopic terms. 

In summary, by combining all the above estimates and by choosing the strength 
of the contact wave 8 CD , the bound on the a priori estimate \ an( i the Knudsen 
number e to be suitably small, we obtain 

M{ - T) + [[ £ \\d a MA)\\ 2 + \\^{u? y \u^){cpX)\\ 2 ]dT 

//«^(.,»,f)^<c E ». 



+ E 

|«|=2 

With the energy estimate, we complete the proof of Theorem 4.1. 
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